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I. INTRODUCTION 

Quantum entanglement means inseparability of joint 
wave function of two or more distant objects into a prod- 
uct of wave functions of individual objects - even in the 
absence of any mutual interaction or communication be- 
tween them. This epitomizes the underlying nonlocal 
character of quantum world. One of the consequences of 
this nonlocal realism is that a single local measurement 
can not reveal the complete state of an entangled sys- 
tem, since the process of measurement itself forces the 
wave function to "collapse" into one its measured (eigen- 
state) state in a probabilistic sense. Thus, measurement 
process can post-selectively play a role in creation and 
manipulation of entanglement, and this is the essence of 
what is called "projective measurement" . 

Efficient generation of entanglement in many-particle 
systems and its robust transmission and transfer to 
other systems is important for quantum information pro- 
cesses. Based on atom-photon interaction and the ex- 
change of photons between the qubits, entanglement in 
distant atomic states 0, d, [|[ and also between pho- 
tons [E HI have been experimentally demonstrated. 
There is another way of entangling two remote systems 
without requiring any direct interaction between them: 
This is based on projective measurement. This indi- 
rect method of creating entanglement between distant 
systems can also be applied for many quantum com- 
munication tasks. In a recent experiment, Moehring 
et al. @ have created entanglement between two dis- 
tant trapped ions by coincident detection of two pho- 
tons spontaneously emitted by the two ions. An ear- 
lier experiment has shown interference of light emitted 
by two atoms @, B @ making use of projective mea- 
surements. Thiel et al. [Io| have proposed a scheme 
of entangling several remote atomic qubits and thereby 
creating Dicke state [ll[ of many-atoms by projective 



measurement of photons using multiple photo-detectors. 
Dicke states are particularly important for their robust- 
ness against particle loss EE LL3| and non-local prop- 
erties of entangled multipartite states [E EE EE Hq |. 
There are several other proposals for projecting distant 
non-interactingparticles into entangled states via photo- 
detection EE EE EE HE HE HI- Continuous variables 
like the quadratures of a field mode (which are analogous 
to position and momentum) have also been employed [23j 
in entanglement studies. 

Bose-Einstein condensate (BEC) is a macroscopic 
quantum object where entanglement arises quite natu- 
rally due to two-body interaction. Bogoliubov theory [24j 
of Bosc condensation reveals that in the ground state of 
condensate, two particles with opposite momentum are 
maximally entangled [25[ in momentum variables as in 
EPR state ;26j . This unique feature makes Bose conden- 
sates a good source of entanglement in motional degrees 
of freedom. Furthermore, in a two-component BEC, one 
can generate ent ang lement in hyperfine spin degrees of 
freedom HE HE HE HE HE HE Q . In order to extract 
the intrinsic entanglement of a BEC for useful purpose 
of quantum information processes, it is required to ex- 
cite quasi-particles in momentum modes by stimulated 
Raman scattering or Bragg scattering [34|, |35|. Then a 
scattered atom becomes entangled with Bragg-scattered 
photon [36[ . In fact, Bragg spectroscopy can be used 
as a tool for generating entanglement of different kinds 
in a variety of physical situations. For instance, tripar- 
tite entanglement among two momentum modes of BEC 
and one electromagnetic field mode can be produced [25[ . 
Furthermore, it has been shown that when a common 
laser beam passes through two spatially separated con- 
densates, photon scattered by the first condensate carries 
and transfers quantum information to the second one and 
thereby two condensates become entangled [33- Similar 
experimental scheme has been used to produce and sub- 
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particles and atom-photon entanglement, we first give a 
brief introduction to Bragg scattering in Sec. II in order to 
reveal the essential physical processes involved in Bragg 
scattering. We then formulate our theoretical model of 
projective measurement and discuss its effect in gener- 
ating entanglement in Sec. III. In Sec. IV, we characterize 
entanglement and discuss our results. The paper is con- 
cluded in Sec.V. 



II. BRAGG SCATTERING 



FIG. 1: A scheme for entangling two separate BECs by Bragg 
scattering and photo-detection. A and B are two BECs which 
scatter photons from pump beams into probe ones. Di and 
D2 are two photo-detectors and 'BS' stands for beam splitter. 



sequently measure phase difference between two spatially 
separated condensates [38j . 

Here we propose a new scheme for generation of entan- 
glement between two remote condensates by projective 
measurement on Bragg-scattered photons. Our scheme 
relies on coincident detection of two Bragg-scattered pho- 
tons coming from two remote condensates in a Han- 
bury Brown- Twiss type experimental arrangement as 
schematically illustrated in Fig.l. A and B are two re- 
mote single-component condensates. Bragg scattering of 
pump photons occur independently at the two conden- 
sates, the scattering is stimulated by probe beams. This 
gives rise to the generation of quasi-particles in both the 
condensates in particular momentum modes determined 
by the relative angle between pump and probe light 
beams. The probe beams are assumed to have the same 
mode, that is, same frequency and polarization property. 
The scattered photons in the probe modes are allowed to 
pass through a beam splitter in order to coalesce them 
to lose which path information. Then the two photons 
coming out of the output ports of the beam splitter are 
detected by two photo-detectors D\ and D 2 in a coinci- 
dent way. The joint projective state of the two spatially 
separate condensates is explicitly non-Gaussian. The en- 
tanglement of such non-Gaussian states can not be ascer- 
tained from the correlation in fluctuations of quadrature 
phase variables or number variables [39[ . To prove the ex- 
istence of entanglement in such a non-Gaussian state, we 
use the inequality criterion recently proposed by Agarwal 
and Biswas [39| based on the phase fluctuation of paring 
operators (of the form ab , where a and b are annihila- 
tion operators of the two su b-sy stems) and the concept 
of partial transpose of Peres [iof and Horodecki [4l[ . Vo- 
gel [42j and also Hillery and Zubairy [43[ have recently 
introduced a class of similar inequalities the violation of 
which is sufficient to show the presence of entanglement 
in two bosonic fields. 

The paper is organized in the following way. Since 
Bragg scattering holds the key for generating quasi- 



In a Bose condensate of weakly interacting atomic 
gases, zero momentum (k = 0) state is macroscopically 
occupied. Therefore, atom-atom collision occurs primar- 
ily between zero- and non-zero momentum atoms. In 
stimulated Raman or Bragg scattering, two far-off reso- 
nant laser beams with a small frequency difference is im- 
pinged on a trapped BEC. There are basically two phys- 
ical processes in Bragg scattering. In the first process, 
a photon from the laser beam with higher frequency is 
scattered into a photon of the other laser mode. This 
causes transformation of a zero-momentum atom into an 
atom of momentum q, where q is the difference in pho- 
ton momentum of the two beams. In the second process, 
an atom moving with a momentum — q is scattered into 
a zero-momentum state. Because of bosonic stimulation, 
the scattering of atoms from zero- to q— momentum state 
will be dominant process. There also occur the processes 
which are opposite to above two processes, but these are 
subdued due to phase mismatch. Thus Bragg scattering 
generates quasiparticles |34| . predominantly in two mo- 
mentum side-modes q and — q. Bragg spectroscopy [j 
with coherent or classical light produces coherent states 
of quasiparticles in a BEC. When these quasiparticles 
are projected into particle domain, they form two-mode 
squeezed as well as entangled state [25| in particle num- 
ber variables. 

Two remote condensates A and B are subjected to 
Bragg scattering with pairs of Bragg pulses. The frequen- 
cies and the directions of propagation of the laser beams 
are so chosen such that Bragg resonance (phase match- 
ing) conditions of scattering in both the condensates are 
fulfilled. We assume the pulse with higher frequency has 
higher intensity and hence can act as a pump. The other 
laser beam which acts as stimulant for scattering is of 
much lower intensity and so can be considered as probe 
beam. We treat pump beams classically. Since scat- 
tering at the two condensates occur independently, the 
Hamiltonian is simply the sum of the Hamiltonian H A 
and H B corresponding to condensates A and B respec- 
tively. Let a q and b q represent the annihilation operators 
for particles with momentum q in condensate A and B, 
respectively. Let the corresponding Bogoliubov quasipar- 
ticles be denoted by a q and (3 q , respectively. In terms of 
these quasiparticle operators, the effective Hamiltonian 
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as derived in the Appendix can be written as 



H eff 



foq (XqXq + X-qX-q) - ^jc)cj 



X-qJ 



H.c. 



(1) 



where the superscript J stands for condensate A or B, 
tj (with j — a,b) is photon annihilation operator for 



the probe beam applied to condensate J. The particle 
operators 7r g (= a q ,b q ) are related to the quasi-particle 
operators x g (= a q ,f3 q ) by Bogoliubov's transformation 



with 



and 



iqXq 



V lX-q 



ftiOq + , 



hut 



,1/2 



(2) 



(3) 



(4) 



is the energy of Bogoliubov's quasiparticle. Here hw q = 
h 2 q 2 /(2m) is the kinetic energy of a single atom, fj, = 

%$- is the chemical potential with £ = (8Tm a a s )~ 1 /' 2 
being the healing or coherence length, 5j is the detuning 
between the pump and probe frequencies, rj = ^/Nf q Q, 



where f q — 



and is the two-photon Rabi 



frequency. We assume the Bragg resonance condition 
((5 ~ uj q ). The hamiltonian can be solved exactly in 
Heisenberg picture. The Heisenberg equations of motion 



for a triad of operators X 

ten in a matrix form X — 
matrix 



a a a_ 



can be writ- 



iujqMX, where M is a 3 x 3 




(5) 



where x — x/uoi 3 . Let D be the diagonalizing matrix of 



M. The solutions can be explicitly written as 



X{t) = DE^D^AXO) 



(G) 



where E 



is a 



diagonal matrix 
?A2t), exp 
AiS being the eigenvalues of M matrix. 



E 



diag.[exp(iAir), exp(iA2r), exp(iA3r)] with r = cu^t and 



III. ENTANGLEMENT PRODUCED BY 
PHOTO-DETECTION 

Our proposed scheme is shown in Fig.l. Quasiparti- 
cles are generated in the condensates A and B due to 
stimulated light scattering in a pump-probe type Bragg- 
spectroscopic method. Let c a and c b denote annihilation 



operators for the two probe light beams scattered by con- 
densates A and B, respectively. Using Eq. (J6]), the scat- 
tered light at the output of the two condensates can be 
represented by 

c a (t) = a q (t)a\ + a- q (t)a- q + a c (t)c a (0) (7) 

c b (t) = b q {t)$l + b- q (t)0- q + b c (t)c b (0) (8) 

where a± q , b± q , a c and b c are time-dependent coefficients 
determined by the Eq. The scattered light output 

coming from the two condensates are passed through a 
beam-splitter and finally collected at the two detectors 
at Dl and D2 as shown in Fig.l. Let the reflectivity and 
transmissivity at left side of the beam-splitter are r and t, 
respectively, while those at right side are r' and t' . Then 
the photon annihilation operators at Dl and D2 can be 
expressed as 



C Dl = fc b 



Cd 2 = tc a + r'c b 



(9) 



(10) 



Let the initial state of the total system i.e. two conden- 
sates plus the two probe fields, be represented by 



*0> =| 0,0) AB | a, P) fields 



fll) 



where | 0,0)^3 indicates a product state with both the 
condensates in ground states of quasiparticles, where first 
'0' corresponds to condensate A and the second '0' to 
condensate B. We assume that both the probe fields are 
in coherent states | a, (3) fields where the field amplitudes 
a and (3 correspond to the probes incident at A and B, 
respectively. Measurement of two-photon correlation via 
coincident detection of scattered probe lights at the two 
detectors will project the two-condensate density opera- 
tor into 



PAB = ■N'^¥ic\daCD 2 CD 1 PoC^ Dl C^ D 



(12) 



where Tr^d^s implies tracing over the field states, po =| 
^E'o)( v I / o I and A/" denotes a normalization factor. Now, 
substituting [5] and [TO] into [TO] and using the relations [7] 
and [5] we obtain pab =\ \ where 



| $> = VAf(a,(3\ C D2 C Dl I * > 

= VA7 {r'f | 0a,5 b (1„2,)) +rt \ S A (l q ,2 q ),0 B )) 

+ VJf(r'r + t't) | E A (l,),Efl(l g )). (13) 

The states | Sj(l q ,2 q )) denotes a superposition state of 
ground, one and two g-phonon excited states of conden- 
sate j(= A,B). Similarly, | Ej(0, 1 9 ),) is another super- 
position state of ground and one phonon excited states 
of condensate 'j'. Explicitly, these superposition states 
can be expressed as 



V2a 2 n | 2 q 
0) A 



+ 2a a a c a I L 



2 2 

+ a r a 



| XU) = a q | l q ) A + a c a | 0) A 



(14) 



(15) 
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Now, we have the reciprocity relations 

r*t' + r't* = 0, r*t + r't'* = 



\r'\ = \r\, \t'\ = \t\, bf + kf = l 

Let t — \t\exp(i<j)) and t' = \t'\ exp(i<j>'). Since phase 
changes by n/2 on reflection, we have r = i\r\ exp(i(j)) 
and r' = i\r'\ exp(i (/>'). Using the reciprocity relations 
and considering the field amplitudes a and (3 as real quan- 
tities, we obtain 

|$> = v0v ? |r||t|(exp(2# / )|OA,5fl)+exp(2^)|5 A ,O s )) 
+ Vj7e X p[i(<f> + <t>')](\r\ 2 -\t\ 2 ) |Ea,S b >. (16) 

For a 50:50 beam splitter, we then have 

I *) =^\{\ 0a,S b ) +exp(2iA0) | S A ,0 B )) (17) 

where Ac/> = <j>- <f/ and | $) = exp(-2i<?!> / ) | $). Eq. ([TTJ) 
is manifestly an entangled state of the two condensates 
A and B. This state obtained via two-photon detection 
is different from the Gaussian state of each independent 
Bose condensate. Furthermore, the basis states involved 
in this entanglement describe collective modes or phonon 
modes of the condensates, | 0a, Sb) refers to a joint con- 
densate state in which condensate A is in zero phonon 
and condensate B is in a superposition of zero, one and 
two phonon states. This state results from quantum in- 
terference of two probable processes which are: (1) two 
scattered photons come from condensate B and no pho- 
ton is scattered by condensate A, these two photons are 
then split by the beam splitter and detected at the two 
detectors projecting the joint condensate state into the 
form | 0,4, Sb}- (2) The second process consists of scat- 
tering of two photons by condensate A and no photon by 
condensate B, this results in the joint condensate state 
I Sa, 0b). Since these two processes are probabilistic, the 
quantum interference of these two processes eventually 
gives rise to the resultant state of the form (fTT)) for a 50:50 
beam splitter. When one scattered photon comes from A 
another from B, there is the probability amplitude of rr' 
that both the detectors will detect the only reflected part 
of the two photons and also the probability amplitude of 
tt' that only transmitted part of both the photons will be 
detected. Both these processes project the joint conden- 
sate state into the form | S^(l g ), Es(l g )). The net prob- 
ability amplitude being the sum of these two quantities, 
for a 50:50 beam splitter they cancel each other. This 
explains why there is no component of | E^(l g ), Eb(1 9 )) 
in Eq. dTTD- 

As discussed earlier, photons scattered by a condensate 
share entanglement with phonons or condensate momen- 
tum modes. These scattered photons can act as the car- 
riers of quantum information of phonons. In a BEC of 
weakly interacting atomic gases, phonons are long-lived. 
By coincident detection of two independent photons but 




FIG. 2: One eigenvalue of the partial-transposed density 
matrix p A ' B calculated in the quasiparticle picture is plot- 
ted as a function of dimensionless interaction time cjf t, for 
different values of average probe photon numbers n p = 10 
(solid), rip — 20 (dashed-dotted). The phase difference is 
A(/> — and the phase difference between the two input laser 
fields Oaf) — 0. The effective atom-field coupling constants 
are 77s = t\a = 7.7uig. 



separately entangled with their respective scatterer con- 
densates, we can establish a quantum communication 
channel between two condensates . Thus our proposed 
method of generating entanglement between two remote 
condensates may also find application in quantum cryp- 
tography [44[ and teleportation (45J. Furthermore, this 
two-photon detection scheme may enable a partial Bell 
state analysis 0, HH, 53] • Most of the earlier proposals 
for creating entanglement between two distant ions or 
atomic ensembles via photo-detection are based on the 
electronic excitation and subsequent emission of photons 
that are detected. In such situations, spontaneous decay 
to unwanted electronic states and decoherence can not 
be avoided. In our proposed scheme, since we use far-off 
resonant stimulated Raman-type light scattering, spon- 
taneous emission is negligible. Since the system is a Bose 
condensate of weakly interacting ultracold atoms with 
long coherence time, decoherence is also at a minimum 
level. Moreover, we use coherent light for stimulating 
photon scattering and so our scheme does not require 
any cavity. 



IV. CHARACTERIZING ENTANGLEMENT IN 
STATE (17) 

A. Negativity of the partial transpose of the 
density matrix 

The necessary and sufficient condition for entangle- 
ment in any bipartite system is the negativity of at least 
one of the eigenvalues of the partial transpose [4(| of the 
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density matrix of the system. Let p A ' B denote the density 
matrix of a bipartite system composed of sub-systems A 
and B. Under partial transpose of Peres and Horodecki 
over the sub-system B (or A), let the density matrix be 
represented by p A ' B (or p A ,B ) which can be derived 
by making transpose only on the operators of B (or A) 
. If the wavefunction of the composite system is insep- 
arable, then there will be at least one eigenvalue of the 
transposed matrix which is negative. 

Now, the state in Eq.(17) can be expressed as 



*) =C | 0,0) 



m— 1 



o | TO,0)+ | C 0m | 0,m)] 



(18) 



where | m, n) represents a joint quasiparticle number ba- 
sis with to number of quasiparticles in condensate A and 
n number of quasiparticles in condensate B. From Eqs. 
(14) and (17), we find 



Cn = a 2 a 2 



exp(2iA(j))b 2 c (3 2 



(19) 



Similarly, all other coefficients C„ l; „ can be deduced from 
Eqs. (14) and (17). The density operator is p AB =| 
|. Taking partial transpose on p AB with respect 
to B implies changing the base operators | ij)(mn |— >| 
in)(mj | and the resulting matrix is p AB . We numeri- 
cally find that in all parameter regimes there is one eigen- 
value of p AB which is always negative and hence the 
state (17) is an entangled state. Figure 2 shows that the 
negativity of the eigenvalue is more prominent when the 



probe photon number \a\ 



is lower. 



B. Violation of the entanglement inequalities for 
observables 



Based on the idea of partial transpose, Simon and also 
Duan et al. [48| have independently given an entangle- 
ment criterion. Using variances in quadrature variables, 
an entanglement parameter 48] can be defined as 



,XP 



(A(X A +X B ) 



(A(P A ~P B ) 2 ) . (20) 



where X s =a,b = (l/V2)(S' + 5 t ) and P s = (l/V2i)(S- 
S^) with S being any bosonic operator of the sub-system 
S. According to the criterion of [48| . the condition for 
the occurrence of entanglement is £xp < 1< This condi- 
tion is necessary and sufficient for Gaussian quadrature 
variables only. For non-Gaussian states, this is only suf- 
ficient. There are certain non-Gaussian bipartite states 
which are conspicuously inseparable but do not fulfill the 
condition £xp < 1- In the present context, the state (17) 
is manifestly a non-Gaussian state. This prompts us to 
look for other criteria [H, |42|, HH based on higher order 
moments of observables. Let us now test whether the 
higher order entanglement criterion introduced in Ref. 
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FIG. 3: Left hand side (LHS) minus right hand side (RHS) of 
the entanglement inequality (|29[) calculated in the quasipar- 
ticle picture is plotted as a function of cof t for average probe 
photon numbers n p — 10 (solid) and n p — 20 (dashed) . All 
other parameters are same as in Fig. 2. 



[39j can reveal the entanglement in state (17). To this 
end, let us first discuss what this criterion is. Using quasi- 
particle operators a and (3 (same as 6t q and [3 q defined 
in Eq.(2), the subscript q is omitted for simplicity), we 
construct the following operators 



K„ 



2i 



tfft -aft) 



K M = - (&a + ftp + 1 



(21) 



(22) 



(23) 



which satisfy SU(1,1) algebra. These SU(1,1) operators 
are previously employed for studying higher order squeez- 
ing. The Heisenberg uncertainty relation of these opera- 
tors implies the inequality 



A 



A 



(24) 



Let us now make partial transpose $ <-> j3>. Under this 
partial transpose , the above inequality becomes 



(&ij3 + d/3t) A 



> (Vd + /3/3 t J> 



(25) 
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After some simple algebra, we have 

2 



A 



(a* (3 + aft 



= N 2 + N + M 



(26) 



aJj3-aft 



= N 2 +N - M - 4|(a t /3)| 2 (27) 



where N 2 = 2{6^aftft), N = (a^a + $ft) and 

M = (d t2 /3 2 ) + (a 2 ft 2 ) - (a f /3 + aft} 2 (28) 

Collecting all these terms, the inequality ([25|) can be ex- 
pressed in the form 

(N 2 + N + M)(N 2 + N - M -4|(a t /3)| 2 ) > \N\ 2 (29) 

Violation of this inequality means the occurrence of en- 
tanglement. It is worth mentioning here that all these 
criteria are sufficient for showing entanglement in a bi- 
partite system. 

We now concentrate on the form of the state given in 
Eq. (fl~8|) . In all the joint base states | m,n), quasiparti- 
cle vacuum state (| 0)5) of either condensate S (A , B) 
appears explicitly. This means that the correlation func- 
tion N 2 of the quasiparticle number operators is zero for 
any parameter regime. This reduces the inequality 
to the form 



M 2 - 4(N + M)\(&* jfy\ 2 ) > 



(30) 



Since N 2 = 0, the variance A (&j3 + aft^ equals to 

(N + M) . Because this variance is non-negative, the 
quantity (N + M) must be non-negative. This implies 
that the inequality (f3"0"]) is violated. Thus we have proved 
that the state (17) is entangled. 



C. Stronger entanglement by selection of probe 
phase 

Here we show that the entanglement can be made 
stronger by choosing appropriate phase of the two probe 
beams. In state (|18p. there is a purely Gaussian com- 
ponent corresponding to the basis | 0, 0) which implies 
that both condensates are in their respective quasiparti- 
cle vacuum. Since the criterion f|29[) is devised for ver- 
ifying entanglement in non-Gaussian bipartite state, we 
expect that, on elimination of this purely Gaussian com- 
ponent from (|18[) . violation of the inequality (f2"9")l should 
be much stronger. This component can be eliminated if 
the coefficient Cq given in Eq. (flT))) is made to vanish. 
For simplicity, let us assume that both the condensates 
are identical and their atom-field coupling constants are 
the same. This means a c = b c . Furthermore, we assume 
that A(f> = 0. In such a situation, to make Co vanish im- 
plies a 2 — — ft 1 where a and (3 are the amplitude of the 
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FIG. 4: The effect of Gaussian component of the entangled 
state (17) on the violation of entanglement inequality (I29p 
in the quasiparticle picture is illustrated. The upper curve 
corresponds to the case when the purely Gaussian basis state 
j 00) is present in (17), and the lower curve corresponds to the 
case when this state is eliminated by choosing the appropriate 
phases of the laser fields. For two identical BEC's and same 
atom-field couplings, when the phase difference between the 
two input probe beams is 7r/2, the coefficient of the basis | 00) 
vanishes provided the path-difference between the two output 
probe beams arriving at the detectors is zero. Here n p — 10 
and all other parameters are same as in Fig. 2. 



two input probe fields which are assumed to be in the co- 
herent states. In other words, this means a — i[3 that is 
the two input probes should have a phase difference 9 a /3 
equal to 7r/2. Now, a phase difference of 7r/2 between 
two beams can easily be made if they are derived from 
a single laser through a 50:50 beam splitter with trans- 
mitted and reflected beams being used as the two input 
beams. We plot left hand side (LHS) minus right hand 
side (RHS) of Eq. (f2"9")l as a function of interaction time in 
Figs. 3 and 4 for various parameters. Obviously, negativ- 
ity of (LHS-RHS) implies violation of the inequality ([30]) 
and hence entanglement. For all our numerical illustra- 
tions, we consider two identically prepared BEC's with 
the same atom-field coupling strength. However, it is 
worth pointing out that this assumption is not necessary 
for exploration of entanglement between the condensates, 
this is assumed only for the sake of simplicity. Figure 3 
shows the effect of probe photon number (or intensity) 
on the degree of violation of the inequality. Weaker the 
probe intensity is, stronger is the violation. Figure 4 il- 
lustrates that when the purely Gaussian component is 
eliminated (9 a p = 7r/2) from the state (17), the negativ- 
ity of (LHS-RHS) becomes much stronger implying that 
the criterion (|2T)|) works better for highly non-Gaussian 
entangled states. For identical coupling constants and 
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FIG. 5: Violation of the inequality (|29[) as a function of cou- 
pling constant r\A = r\B = i) for two cases 9 a /3 = (solid) and 
da/3 = (dashed) which correspond to presence and absence 
of the purely Gaussian basis | 00) in the entangled state. Here 



ui q t = 5 and n p = 10. 



detuning between the pump and the probe for the two 
condensates, we find that the degree of violation of the 
inequality saturates at the long interaction time limit. 
However, for any mismatch in coupling constants or de- 
tunings, we have found that (LHS-RHS) tends to zero 
in the long interaction time regime (not shown). The 
effect of atom-field coupling strength on the negativity 
of (LHS-RHS) is illustrated in Fig. 5. We notice that 
in the strong-coupling regime the negativity is stronger. 
Furthermore, in the limit of very strong-coupling the vio- 
lation of the inequality becomes almost insensitive to the 
initial phase-difference d a $ of the two probe beams. This 
indicates that strong-coupling is important for entangle- 
ment between remote condensates. This in turn brings 
in the important role of condensates in quantum infor- 
mation science: Since the effective atom-field coupling 
strength is proportional to the square root of number 
of atoms, strong-coupling regime can easily be attained 
with a BEC rather than a thermal gas. Since thermal 
and phase fluctuations in a BEC is at the minimum level, 
collective atom-field coupling can easily be accomplished 
with a BEC in a cavity with a moderate Q-factor. For 
experimental verification of entanglement, one needs to 
measure the various variances (appearing in the inequal- 
ity) in the phonon or collective excitation modes of BEC. 



V. CONCLUSION 

In conclusion, we have shown that two remote inde- 
pendent condensates can be made entangled in collective 



excitations of BEC such as quasiparticle or phonon vari- 
ables by projective measurement on two photons Bragg- 
scattered by the two condensates. The generated en- 
tangled state is explicitly non-Gaussian and the existing 
criterion for entanglement in Gaussian variables is insuf- 
ficient in revealing entanglement in this state. This has 
prompted us to test another criterion [39[ introduced spe- 
cially for testing entanglement in non-Gaussian states. 
We have shown that with this criterion the projective 
state of the two condensates is entangled and the entan- 
glement can be made stronger by choosing probe phase in 
such a way so that the non-Gaussian nature of the entan- 
gled state becomes more prominent. The entanglement 
is shown to be stronger when the probe beams are weaker 
and the atom-field coupling is higher. Since atoms behave 
collectively in BEC-field interaction and also because of 
yN (N is the atom number) scaling of coupling strength, 
BEC can play an important role in generation and ma- 
nipulation of entanglement. Since the entanglement in 
the present context arises in the higher order fluctua- 
tions of the paring operators, experimental detection of 
this entanglement requires new techniques that will en- 
able to measure the pairing fluctuations of phonons or 
quasiparticles. Using Bragg spectroscopy phonons in a 
BEC have been detected by Ketterle's group. But how 
to detect phonon-phonon correlation which is crucially 
required for exploration of higher order entanglement is 
presently unknown. Perhaps by knocking out pairs of 
atoms from the two trapped condensates and measur- 
ing their fluctuations, one can detect this entanglement. 
Our proposed scheme can be easily generalized for mul- 
tiple condensates. Furthermore, it may be interesting to 
generate entanglement both in spin and motional degrees 
of freedom using spinor condensates and polarized lights 
in our scheme with the possibility of interesting inter- 
play of entanglement in spin and center-of-mass degrees 
of freedom. 
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APPENDIX 

Here we derive the effective hamiltonian for Bragg scat- 
tering in a condensate. The total hamiltonian of a con- 
densate interacting with two single-mode light fields is 
H = Ha+Hf+Haf, where Hp = huiic^c^ +foui2C] l2 C] si2 
corresponds to the two fields described by the operators 
Ck x and Ck 2 with photon momentum ki and k2 and the 
frequencies lo\ and 0J2, respectively. We assume uj\ > lu2- 
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The free part of the atomic hamiltonian 



7T k 7T k 



k 



2mV 



X ^k 3 ^k 4 ^ k ^ k 6' 5k 3+k4:k5+k 6 (31) 

k3,k 4 ,k 5 ,k6 

governs the dynamics of a weakly interacting atomic con- 
densate and 

H A F = ^4 2 Ck! X] (^+k^k + TT-q+kTTk 



H.C. 



(32) 



describes atom- field interaction. Here TTk(Tt k ) is the 
annihilation(creation) operator of an atom with mo- 

mentum k and frequency ujk — q = ki — k2, 

Q = (_Ei.di 3 )(i?2-rf32)/(^ 2 A) is the two-photon Rabi fre- 
quency, where £?i(2) are the field amplitudes, the dij is the 
electronic transition dipole moment between the states 
\i) and \j) of an atom and A is the detuning of the first 
laser field (with frequency from the transition fre- 
quency between the electronic ground (|1)) and excited 



(| 3)) levels of the atom. For a single-component conden- 
sate the electronic ground states |1) and |2) are the same. 
Here a s is the s-wave scattering length of the atoms and 
V is the volume of the condensate. 

Using Bogoliubov's prescription ttq , ttq — > \/No, and 
keeping the number density hq = Nq/V fixed in the 
thermodynamic limit, one can transform the hamiltonian 
Ha into a quadratic form. Further, applying Bogoli- 
ubov's transformation it is possible to diagonalize Ha 
and rewrite the entire hamiltonian in terms of Bogoli- 
ubov's quasi-particle operators Xk- Considering the con- 
densate ground state energy as the zero of the energy 
scale, and treating the laser light with higher frequency 
(uji) classically, the effective hamiltonian can be written 



H, 



eff 



XqXq + X-qX-q) ~ ^4 2 ^ 



^ £ k 2 (Xq + X-q) +H.C. 



(33) 



where S = uj\ — u>2 and r\ = y/~Nf q Cl; where f q = u q — 
v q . In writing the above equation, we have retained only 
two dominant momentum side-modes of the condensate 
under the assumption of Bragg resonance (5 ~ uo q ). 
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